Emergent Critical Phase and Ricci Flow in a 2D Frustrated Heisenberg Model 
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We introduce a two-dimensional frustrated Heisenberg antiferromagnet on interpenetrating hon- 
eycomb and triangular lattices. Classically the two sublattices decouple, and "order from disorder" 
drives them into a coplanar state. Applying Friedan's geometric approach to nonlinear sigma mod- 
els, we obtain the scaling of the spin-stiffnesses governed by the Ricci flow of a 4D metric tensor. 
At low temperatures, the relative phase between the spins on the two sublattices is described by a 
six-state clock model with an emergent critical phase. 

PACS numbers: 75.10.-b, 75.10.Jm 



A remarkable discovery of recent years is that frus- 
trated two dimensional Heisenberg models can evade the 
Hohenberg-Mermin- Wagner theorem [1, 2] via the devel- 
opment of long-range discrete order driven by short-range 
thermal spin fluctuations: such discrete long-range order 
develops despite the persistence of a finite spin correla- 
tion length, leading to a finite temperature Ising (Z2) or 
Z3 Potts phase transition [3-9] . This phenomenon is well- 
established in the J\-Ji Heisenberg model on the square 
lattice and has recently been realized in iron-based su- 
perconductors [10]. An interesting question motivated by 
this discovery is whether it can be generalized to higher 
(p > 5) order. If one can show, in addition, that these 
emergent discrete degrees of freedom are described by a 
p-state clock model [11, 12], the unique situation arises 
that a Heisenberg spin system exhibits two Berezinskii- 
Kosterlitz-Thouless (BKT) transitions which bracket a 
critical phase. In a system of discrete Ising spins, such 
a scenario was reported to occur on the triangular lat- 
tice [13, 14]. 

In this Letter, we introduce such a Heisenberg model 
defined on interpenetrating honeycomb and triangular 
lattices (Fig. 1(a)) with nearest- neighbor antiferromag- 
netic coupling. This model may be realized with cold 
spinful atoms in optical lattices, where it arises naturally 
in the limit of large on-site interactions [15-18]. An- 
other promising experimental route is to employ STM 
techniques for nano-fabrication and spin-resolved read- 
out of stacked triangular and honeycomb monolayers of 
magnetic atoms like Cr or Co [19-22]. For classical 
spins the two sublattices are decoupled giving rise to an 
SO(3)xO(3)/0(2) order parameter. "Order from disor- 
der" [23, 24] drives the two sublattices into a coplanar 
spin configuration [25] with an SO(3)xU(l) order pa- 
rameter and a six-fold in-plane potential. In the coplanar 
state we explicitly show that the U(l) degrees of freedom 
decouple to form an emergent Zq clock model with an in- 
termediate power-law phase. This non-trivial decoupling 



of the U(l) phase is essential for the critical phase to 
occur. 

A novel aspect of our work is that we apply Friedan's 
coordinate-independent approach to nonlinear sigma 
models [26] to the scaling of the spin-stiffness. In this 
approach the configurations of the 2D spin system cor- 
respond to a worldsheet of a string evolving in four- 
dimensions, where the metric is determined by the com- 
ponents of the antiferromagnetic stiffness and its renor- 
malization corresponds to a Ricci flow of the metric ten- 
sor. We also note the decoupling of our U(l) phase can 
be viewed as a toy model for the compactification of a 
four-dimensional string-theory. 

Specifically we study the antiferromagnetic Heisenberg 
model on a decorated 2D triangular lattice (cf. Fig. 1(a)); 




FIG. 1. (color online) (a) Heisenberg model on the "windmill" 
lattice, (b) Definition of angles a and ft describing the relative 
orientation of magnetic order on triangular and honeycomb 
lattice, (c-d) Angle dependent free energy correction SF from 
thermal and quantum spin fluctuations for parameters J^h = 
Jtt = 1, Jth = 0.4, T = 1. Panel (d) is for fixed /3 = tt/2. 
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the associated Bravais lattice has three basis sites per 
unit cell at positions b t = ao(0, 2/v / 3), bA = (0,0) and 
&b = ao(0, l/VS) where indices A, B label the two hon- 
eycomb sites. We set the lattice constant a = 1. The 
Hamiltonian is H = H u + Hab + UtA + HtB with 
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where S a (rj ) denote spin operators at Bravais lattice site 
j and basis site a £ {t, A, B}. The vectors {S a b} point 
between nearest-neighbors of sublatticcs a, b. We assume 
in the following that the spin exchange couplings within 
the same sublattice are larger than the inter-sublattice 
coupling J th < J t t,Jhh, where J t h = JtA = JtB and 
J hh = Jab- For decoupled lattices J t h — 0, the clas- 
sical ground state on the bipartite honeycomb lattice is 
the usual Neel state, while spins on the triangular lat- 
tice arrange in a 120°configuration [27]. Although the 
exchange fields between the two sublattices exactly can- 
cel for this configuration even for J t h > 0, quantum and 
thermal fluctuations depend on the relative orientation of 
the magnetization on the two sublattices. The uniaxial 
magnetic order on the honeycomb lattice is described by 
a normal vector n(x), which points along the magnetiza- 
tion on sublattice A. The biaxial order on the triangular 
lattice is characterized by a triad of orthonormal vectors 
{tj(x)} with j = 1,2,3. Equivalently, it may be ex- 
pressed by an orthogonal matrix t = (tx,t%,t^) G 50(3). 
We take the vectors £1.2 to span the plane of the mag- 
netization on the triangular lattice. The relative order 
between the two sublattices is thus determined by two 
angles a and /3, that are defined in Fig. 1(b). 

Symmetry considerations dictate the form of the long- 
wavelength action which takes the form of a nonlinear 
sigma model (NLSM) 
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The action contains the usual gradient terms of the 
0(3)/0(2) and the 50(3) NLSM for the order param- 
eter on the honeycomb and triangular lattice. The bare 
spin stiffnesses K,Kj can be derived in a 1/5-expansion 
and read K = 2J hh S 2 /T, K 1 = K 2 = V3J tt S 2 /AT and 
K 3 = [28 30]. In addition, the action in Eq. (2) con- 
tains two potential terms, generated by short-wavelength 
spin fluctuations ("order from disorder") [23, 24] 

S c = \ J d 2 x (7 cos 2 (3 + A sin 6 (3 sin 2 (3a)) . (3) 

A positive 7 > favors coplanarity, whereas 7 < in- 
duces n to be perpendicular to the plane of the triangular 
magnetization. The six-fold anisotropy term A is relevant 
only for 7 > 0. 
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FIG. 2. (a) Schematic phase diagram, (b) Coplanar RG flow 
of the variables l' a (blue, increasing), / = (Iihh) 1/:i (greed 
dashed), (I2 — Ii)/I (red), and r (pink dotted). Curves are 
normalized to initial values at Z 7 . Upper panel is for J t t 2> Jhh 
with Jhh = 1, Jtt = 5, Jth ~ 0.4, T = 0.6, and initial values 
/ = 5.3, {h -h)/I = 0.27, r = 0.82, I' a = 1.2. Decoupling 
is due to (Ji — — > 0. Lower panel is for Jhh 2> Jtt 
with Jhh = 5, Jtt = 1, Jth = 0.4, T = 0.5, and initial values 
I = 4.5, (I 2 - h)/I = 2.1, r = 0.11, I' a = 1.1. Decoupling is 
due to r — > 0. 



Heuristically, we expect 7 > 0, favoring coplanarity: 
spins on the honeycomb lattice can minimize their en- 
ergy by aligning themselves perpendicular to the fluc- 
tuation Weiss field from the triangular lattice [24]. To 
confirm this reasoning, we have performed a Holstein- 
Primakov spin wave analysis of Eq. (1). Our results for 
the fluctuation correction to the free energy for arbitrary 
angles a and /? between the two sublattices are given in 
Fig. l(c-d) and show that 7 > 0. For small Jth we find 
7 = (Jth/ J) 2 A 1 (Jtt/ Jhh, J /T) is the dominant term 
in the potential, while A = (J t h/J) 6 Ax (Jtt/ Jhh, J /T) , 
where J = \J Jtt Jhh and the ^4 7 ,a are functions that de- 
pend weakly on J u /Jhh- 

As temperature is reduced, the two sublattices enter a 
coplanar regime. The temperature scale for this crossover 
is easily determined from standard scaling arguments and 
yields the coplanar crossover temperature (see Fig. 2(a)) 
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in case where Jhh < Jtt- In the opposite regime Jtt < Jhh 
we obtain an implicit expression for T cp that also ap- 
proaches zero only logarithmically for 7 — >■ 0. The 
crossover temperature in Eq. (4) follows from the known 
flow equation — — l/2n for the stiffness K with run- 
ning cutoff A(l) = aQ 1 e~ L and the flow of the coplanar 
potential amplitude j(l) = 7exp(2Z) that is determined 
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by its engineering dimension. While the spin stiffnesses 
are reduced at longer length-scales, the potential term 
grows, and scaling stops when 7(Z 7 ) = 1, which defines 
a length-scale a 7 = a e l ~i ~ a {J hh / J th ) 2 . The coplanar 
crossover takes place when this length-scale is compa- 
rable to the shorter of the magnetic correlation lengths 
on the two sublattices. From the known flow equation 
of the 0(3) /0(2) and the 50(3) NLSM further follows 
that the stiffnesses of the triangular lattice approach an 
isotropic fixed point [31]. The six- fold symmetric poten- 
tial oc A flows to larger values, yet due to A <C 7 holds 
that A(Z 7 ) ~ 0(Jf h ) < 1. 

Once the two sublattices are coplanar, their dynam- 
ics are intimately connected. To describe this regime we 
impose a hard-core constraint: n _L f 3 , i.e., (3 = x/2. 
It is now convenient to introduce a second triad /ii,2,3j 
defining an 50(3) matrix h = ( y hx,h2,h 3 j that de- 
scribes the magnetic order on the honeycomb lattice 
with h\ = n. The coplanar constraint is expressed as 
t = hU where U = exp(iaT 3 ) determines the relative in- 
plane orientation of the two sublattices, defined by the 
angle a. We describe h in terms of three Euler an- 
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Here, the r a satisfy the 



SU(2) algebra [r Q , r^] = ie a bcT c and take the adjoint form 
(T a )bc — i^bac- The coplanar system is thus determined 
by an 50(3) x U(l) order parameter, defined by three 
Euler angles and a single relative phase a. 

To analyze this coupled problem we write the action 
in the form 5 = Sx + 5 C , where 

S * = \f d 2 xg ij [X(x)]d^X i (x)d fl X^x) (5) 

with coordinates X = (<p, (9, ip, a) and stiffness tensor 
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with Sxi — sinX 1 and c^i = cosX J . In our system we 
find h = K x +K 3 , h = Ki + K 3 + K, I 3 = 2K X + K, 
which are set by the stiffnesses of the two sublattices at 
I = Z 7 . The £7(1) degree of freedom a has an initial 
stiffness I a = 2K\{1 1 ) and is coupled to the non-Abelian 
50(3) sector by the term JC = f(cg,0, 1) in the four- 
dimensional metric, where k = 4Kx(l~). The six- fold 
potential 5 c (/3 = |) = |A/ d 2 x sin 2 (3a) is a small but 
relevant perturbation to Sx- At length-scales where A is 
small, the anisotropy 5 C and the gradient term Sx (5) is 
the action of a classical string in a four dimensional space 
with coordinates X (x) at the two dimensional world- 
sheet point x, with metric tensor ^ [X]. Under coordi- 
nate transformations Xi — > X[, Sx in Eq. (5) is invari- 
ant, with transformed metric g' lm = gij ox' m • Like 



Einstein's theory of gravity, this covariance tells us that 
the long-wavelength action Sx is co-ordinate indepen- 
dent and only depends on the geometric aspects of the 
mapping X (x) of the wordsheet to the compact four- 
dimensional space of the coordinate X. The rcnormaliza- 
tion group (RG) flow of the metric tensor must also be co- 
variant under co-ordinate transformations, and following 
the geometric formulation of the NLSM by Friedan [26], 
to two loop order takes the form 
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where R 



iklr, 



is the Riemann curvature tensor and i?„ = 



R k ikj is the Ricci tensor [32]. This expression defines a 
generalized Ricci flow [33]. The Riemann tensor is de- 
termined by the Christoffel symbols T l - k = \g ll {gji,k + 

9ki,j - 9jk,i) as R k uj = rf j4 - rfcj + r^rg - r^.rg. 

The flow equations of our five coupling constants Ij, I a 
and k follow from Eq. (7). 

A key insight into the low energy phase diagram is ob- 
tained by noting the coupling term K, can be eliminated 
via a coordinate transformation tp tp' = tp + ra with 
r = k/2I 3 . This yields a metric g in Eq. (6) with K = 0, 
I a I' a = I a — k 2 /4/ 3 yet with g so W that depends on 
the U (1) phase a via the above shift of the Euler angle tp. 
This gauge transformation to the appropriate center of 
mass coordinates allows for clear criteria when the U(l) 
sector of the theory decouples from the 50(3) sector: if 
either \h - I 2 \ < \flxh or r « 1 it follows that g so ^ 
becomes independent of a and the 17(1) phase decouples 
from the dynamics of the non-collinear magnetic degrees 
of freedom. The first criterion follows from the fact that 
gSO(3) j g independent of tp if I\ = 1%, while the second 
criterion implies that the shift in ip is negligible. From 
Eq. (7) follows after a lengthy but straightforward calcu- 
lation that 71,2,3 flow t° an isotropic fixed point, while 
the dimensionless variable r follows the flow equation (for 
simplicity we only list the one loop result, the two loop 
correction does not change our conclusions): 
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Thus, if the initial anisotropy \I% — J2I — K is weak, 
which happens for Jhh *C Jtt, the coupling r does not 
change much. The 50(3) sector, however, quickly be- 
comes isotropic in the 1-2-plane leading to a decoupling 
of the Z7(l) phase. On the other hand, in the limit of 
strong anisotropy for 3> Jtt, where \I\ — /2I is not 
small, we find that r vanishes rapidly. In both cases fol- 
lows that the phase angle a emerges as an independent 
degree of freedom. The /3-function for the reduced phase 
stiffness I' a = I a — k 2 /4/ 3 follows from Eq. (7) as 
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and does, as expected, approach zero once either of the 
two decoupling conditions are fulfilled. Thus, perturba- 
tively no renormalization of the stiffness I' a takes place. 
In Fig. 2(b) we present the coplanar renormalization 
group flow for two different sets of parameters corre- 
sponding to weak and strong initial anisotropy. An in- 
teresting aspect of the decoupling follows from the Ricci 
scalar R = g^R^: 
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is the Ricci 



2/1/2/3 i, 

scalar of the SO (3) sector. Once the decoupling takes 
place, P a — > and the U (I) sector becomes flat. On 
the other hand R -> R so <- 3 ) grows under renormaliza- 
tion since the stiffnesses Ij decrease. Thus, we arrive 
at a flat one dimensional sector weakly coupled to an 
three-dimensional manifold with large curvature. This 
"curling-up" and asymptotic decoupling of a subspace 
may serve as a toy model for compactification. 

Since the decoupling emerges rapidly in both limits 
Jhh *C Jtt and Jhh ^ Jtt, we find that A, whose flow 
is governed by ^A = (2 — 9/nI' a )\, is still small at the 
decoupling lengthscale. The resulting low-energy theory 
corresponds to S = Sso(3) + Sze with 



S ^ = \J d 2 x[(I' a (d^a) 2 + \ S m 2 (3a)] 



(11) 



This is the well-known six state clock model that ex- 
hibits two consecutive BKT transitions [11]: one at 
Tg KT that separates a high temperature disordered phase 
from a low temperature critical phase, where correla- 
tions (exp[z(a(x) — a(x'))]) decay with a power-law in 
\x — x'\; a second at T^ KT where the symmetry is 
spontaneously broken, leading to true long-range order 
with a = mr/3 (n <E {1,...,6}). It is crucial that the 
decoupling of the U(l) phase occurs first, otherwise the 
50(3) sector would screen the long-range interactions be- 
tween topological defects - vortices at T^ KT or domain 
walls at T^~ KT that are responsible for the BKT transi- 
tions and the intermediate critical phase. 

Following the RG program of the BKT problem for 
Eq. (11) we need to take into account that the size of 
the vortex is now given by the coplanar lengthscale a 7 3> 
ao [34, 35]. We determine the vortex unbinding transition 
temperature T^ KT implicitly via 
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with fugacity y — e 5c a^/a§ and core action S c ~ n{l + 
mm(K,K t )}. From Eq. (12) we predict that T> KT < 
T cp , i.e., the BKT transition is only numerically smaller 
than the coplanar crossover temperature. The system 
enters the critical phase soon after it becomes coplanar. 



Similarly, it follows from Ref. [11] that T^ KT and Tg KT 
are of the same order of magnitude. The resulting phase 
diagram is shown in Fig. 2(a). 

In summary we have presented a 2D Heisenberg model 
on a decorated triangular lattice where short wavelength 
thermal fluctuations select long-range Z@ order. This is 
preceded in temperature by an emergent critical phase 
that is framed by two BKT transitions. We have written 
the action of this model as a classical 4D string theory 
where the spin stiffness is determined by the metric ten- 
sor of the manifold; the scaling equations are then ex- 
tracted as components of the resulting Ricci flow. Wc 
note that the decoupling of the £7(1) degree of freedom 
corresponds to a dimensional reduction of the analogous 
string theory and thus to a toy model of compactification. 
Finally we note that the emergence of massless modes in 
collective mode massive theories could have interesting 
implications for two-dimensional field theories. 
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Note added.- After obtaining these results we learned 
of two recent studies: one on a Kitaev-Heisenberg model, 
where an emergent Zg-symmetry results from a concep- 
tually different mechanism [36] , a second on itinerant sys- 
tems where an emergent Z4 Potts model appears [37]. 
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